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ABSTRACT 

Neutralinos may be captured in the Sun and annihilate therein producing 
high-energy neutrinos. Present limits on the flux of such neutrinos from un- 
derground detectors such as IMB and Kamiokande II may be used to rule out 
certain supersymmetric dark-matter candidates, while in many other supersym- 
metric models the rates axe large enough that if neutralinos do reside in the 
galactic halo, observation of a neutrino signal may be possible in the near fu- 
ture. Neutralinos that are either nearly pure Higgsino or a Higgsino/gaugino 
combination are generally captured in the Sun by coherent scattering off nuclei 
via exchange of the lightest Higgs boson. If the squark mass is not much greater 
than the neutralino mass then capture of neutralinos that are primarily gaugino 
occurs predominantly by spin- dependent scattering off hydrogen in the Sun. The 
neutrino signal from annihilation of WIMPs with masses in the range 80-1000 
GeV in the Sun should generally be stronger than that from WIMP annihilation 
in the Earth, and detection rates for mixed-state neutralinos are generally higher 
than those for Higgsinos or gauginos. 
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I. INTRODUCTION 
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The idea that stable weakly-interacting massive particles (WIMPs) make up 

the bulk of the dark matter in the Universe and in the galactic halo has been 

the focus of much theoretical and experimental research recently. 1 Now that the 

original WIMP, the Dirac neutrino, has been ruled out, the neutralino — a 

linear combination of the supersymmetric partners of the photon, Z°, and Higgs 

bosons — has become the preferred thermal relic. Although the original treatises 

considered only neutralinos lighter than the W^, 4 ’ 5 heavy neutralinos — those 

more massive than the W — may also be suitable dark-matter candidates. ' Al- 
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though “extremely” massive neutralinos are not favored theoretically, neutrali- 
nos in the 100-GeV range may still solve the naturalness problem and become 
increasingly attractive as unsuccessful accelerator searches push the mass scale 
for supersymmetry upward. 

Since many neutralinos are not yet accessible in accelerators and are such 
compelling dark-matter candidates, a variety of complementary experiments to 
detect neutralinos in our galactic halo are currently being pursued. Some seek to 
observe neutralinos by detecting the energy deposited in am ultra-low background 
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detector when a neutralino elastically scatters off of a nucleus therein. Alterna- 
tively, neutralino dark matter in the galactic halo may be indirectly detected by 
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its annihilation products. A continuum spectrum of cosmic-ray antiprotons, 7 
rays, and positrons, are produced in the cascade resulting from the annihi- 
lation products of the neutralinos; however, astrophysical uncertainties involving 
the propagation of cosmic rays from conventional sources are so great that it 
seems unlikely that WIMP-induced continuum cosmic rays could ever be distin- 
guished from those from standard sources. Some authors have boldly suggested 
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that annihilation of WIMPs in the galactic halo could produce either 7-ray 
or positron 14,15 line radiation which could be readily distinguished from back- 
ground. While such a signal would provide unambiguous evidence for particle 
dark matter, because of astrophysical uncertainties an observable signal of this 
kind is not guaranteed even if suitable WTMPs do reside in the galactic halo. 

In this paper we address the possibility of indirect detection of heavy neutrali- 
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nos by observation of yet another annihilation product: high-energy neutrinos. 
WIMPs in the galactic halo will be captured in the body of the Sun or the 
Earth 16,17,18 and annihilate therein producing high-energy neutrinos that may 
be observable in underground neutrino detectors. This method of detection has 
several advantages over cosmic-ray signatures: First of all, whereas cosmic rays 
are expected to be isotropically distributed, the neutrino signal comes from a 
fixed direction and is therefore much more easily distinguished from background. 
The number density of neutralinos in the halo is inversely proportional to 
the neutralino mass and, as we shall see, the annihilation rate in the Sun is 
oc n* while the annihilation rate in the halo is oc n* 2 , making the neutrino 
signed favored for higher neutredino masses. In addition, the uncertainties in 
the predicted rates for neutrino events are smaller than those in the predicted 
cosmic-ray fluxes (roughly factors of about two for neutrino events and orders of 
magnitude for cosmic-ray fluxes). Basically this is because the local halo density 
is known better than the dark-matter distribution throughout the galaxy, and 
propagation of neutrinos through the Sun is more easily modeled than cosmic-ray 
propagation through the galaxy. It should also be noted that neutrino searches 
and cosmic-ray searches Eire mutually complementary: For example, the neutrali- 
nos that may be discovered through distinctive cosmic-ray positron signatures are 
primarily Higgsinos , 15 whereas neutrino signals are strongest for neutralinos that 
are a mixed Higgsino/gaugino state. 

Unlike Dirac neutrinos which annihilate directly into light (».e., v t , and 
u T ) neutrinos, neutralinos are Majorana particles and therefore do not produce 
prompt neutrinos; the neutrinos from neutralino-neutralino annihilations come 
from the decays of the annihilation products, so the neutrino spectrum is con- 
siderably softer. Detailed neutrino spectra from energetic quarks and leptons 
V 19 

injected into the core of the Sun were calculated by Ritz and Seckel (RS). The 
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analysis for light neutralinos was originally carried out by Giudice and Roulet 
who considered only annihilation into fermion- antifermion pairs and more com- 
pletely by Gelmini, Gondolo, and Roulet who considered annihilation into 
pairs of Higgs bosons as well. Here we extend this work to heavy neutralinos by 
considering the effect of the gauge-boson, Higgs-boson, and top-quark annihila- 
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tion channels which open up for heavy neutralinos. We also consider the effect 
of interactions of the annihilation products and resulting high-energy neutrinos 
in the Sun which become important at higher energies. 

First let us briefly review the minimal supersymmetric extension of the Stan- 
dard Model (MSSM) and the properties of the neutralino. For more details we 

6 

refer the reader to Ref. 3 and Griest, Kamionkowski, and Turner (GKT) whose 
notation we use throughout. There are actually four neutralinos, and the lightest 
(the nth) is assumed to be the lightest supersymmetric particle (LSP) and stable 
and is denoted as the neutralino, 

X = ZnlB + Z n2 W 3 + Z„3#l + Z n4 ff 2 , (1) 

where ( Z)ij is a real orthogonal matrix that diagonalizes the neutralino mass ma- 
trix [Eq. (C38) in Ref. 3] and depends only on the gaugino mass parameter M, 
Higgsino mass parameter n, and the ratio of Higgs vaccuum expectation values 
tan / 3 . In Fig. 1 we plot neutralino mass contours (broken curves) and contours of 
Z 2 i + Z 2 2 (solid curves), the gaugino fraction, for tan (3 = 2 (plots for other values 
of tan/? are similar). As noted originally by Olive and Srednicki in much of 
parameter space where the neutralino is heavier than the W , the gaugino fraction 
is greater than 0.99 and the neutralino is almost pure B-ino. In much of param- 
eter space, the gaugino fraction is less than 0.01 and the neutralino is almost 
pure Higgsino. Near the 0.5 gaugino fraction curve, a curve that asymptotes to 
fj, = |M tan 2 Bw at high neutralino mass, the neutralino is a mixed state, half 
gaugino and half Higgsino. 

In the MSSM there axe three neutral Higgs bosons. The mass of the lightest, 
#2 — which must be less than mg cos 2/3 (provided the top quark is not unusually 
heavy; see Ref. 23) — and tan 0 determine the masses of the other two, H ® — 
which must be heavier than the Z — and H 3 — whose mass falls between m # 0 and 
m fjo. There are also charged Higgs bosons which are always heavier than the 
W and two charginos, linear combinations of the supersymmetric partners of the 
W and charged Higgs bosons. The masses of the superpartners of the quarks and 
leptons, which we will collectively refer to as squarks, axe all undetermined, but 
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for simplicity we give them all the same mass Mg which, assuming the neutralino 
is the LSP, is greater than m*. 

3 

Although the MSSM has many undetermined parameters (tan 3, M, /x, 
rriffo, Mq, and the top-quark mass m<) the parameters are not entirely uncon- 
strained, and by studying several “comers” of parameter space we can get an 

understanding of the dependence of detection rates on the different parameters 
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of the model. Although mj is constrained only to be greater them 80 GeV (from 
unsuccessful accelerator searches) and less than about 200 GeV (from limits 
on radiative corrections to sin 2 $w), we will assume m< = 120 GeV throughout; 
as we will discuss later, varying the top-quark mass should have little effect on 
our results. Recent searches for neutral Higgs bosons at LEP have constrained 
regions of m //o-tan {3 space. In addition, we will only consider tan/3 > 1, since 
radiative corrections drive tan/3 to values greater than one when m* mj, and 
tan/3 < mj/mj ~ 25, required for electroweak symmetry breaking in many su- 
pergravity models. To see the range of possible capture and detection rates 
due to the range of all possible values for the squark mass we will present results 
assuming the squark mass is infinite and then show results assuming the squark 
mass is slightly heavier than the neutralino mass. 

Although determination of the event rate is relatively straightforward, it is 
quite lengthy and depends on a variety of input physics such as solar physics, 
neutrino physics, hadronization of quarks, underground detectors, and, of course, 
the interactions of neutralinos with ordinary matter. The flux of high-energy 
neutrinos of type i (e.g., i = etc.) from neutralino annihilation in the Sim 

is simply 



The quantity is the rate of neutralino-neutralino annihilations in the Sun, 
and R is simply the distance of the Earth from the Sim. Neutralinos from the 
galactic halo are accreted onto the Sun and their number in the Sun is depleted 
by annihilation. In most cases of interest these two processes come to equilibrium 
on a time scale much shorter than the solar age in which case = C/2 where 
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C is the rate for capture of neutralinos from the halo. As one might imagine, the 
capture rate is basically determined by the flux of neutralinos incident on the Sim 
and a probability for capture which in turn depends on kinematic factors and the 
cross sections for elastic scattering of the neutralino off of the elements in the Sun. 
The sum is over all annihilation channels F (e.g., pairs of gauge or Higgs bosons 
or fermion-antifermion pairs), Bp is the annihilation branch for channel F, and 
( dN/dE)pi is the differential energy flux of neutrino type i at the surface of the 
Sun expected from injection of the particles in channel F in the core of the Sun. 
The flux ( dN/dE)fi is a function of the energy of the neutrino and of the energy 
of the injected particles. Determination of these fluxes is quite complicated as it 
involves hadronization of the annihilation products, interaction of the particles 
in the resulting cascade with the solar medium and the subsequent interaction of 
high-energy neutrinos with the solar medium as they propagate from the core to 
the surface of the Sun. Neutralinos may also be captured in the Earth; however, 
for a number of reasons which we will discuss below, the rates for neutrino events 
from neutralino annihilation in the Earth will generally be smaller than those 
from the Sun if the neutralino is heavy. 

The experimental signature on which we will eventually focus will be the 
number of upward-moving muons induced by high-energy neutrinos from the 
Sun that are observed in underground detectors. Given the fluxes ( d<j>/dE)i the 
final result for the rate (per unit detector area) for neutrino-induced upward 
moving muons may be written simply as 

^detector = ^ A j E 2 dE, (3) 

where the sum is over which produce muons, and £7,, which produce an- 
timuons. Since the cross section for the neutrino to produce a muon in the rock 
below the detector is proportioned to the neutrino energy E and the range of the 
muon is roughly proportional to its energy, the probability a neutrino of energy 
E produces a muon which traverses the detector is F 2 times a constant Dp, hence 
the integral in Eq. (3). Neutrinos may also be detected by contained events in 
which a charged lepton is produced within the detector, but because this process 
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is proportional only to the neutrino energy E (as opposed to E 2 for throughgoing 
events), the throughgoing muons should provide a more promising signature for 
heavy neutralinos. 

In the next Section we discuss the rate of neutralino-neutralino annihila- 
tion in the Sun and in the Earth. The annihilation rate is proportional to the 
square of the number of neutralinos in the Sun or Earth, and this number is in- 
creased by capture of neutralinos from the halo while neutralinos are depleted by 
annihilation. Capture occurs by elastic scattering of neutralinos in the galactic 
halo off of nuclei in the Sun. . We show the regions of parameter space in which 
capture occurs predominantly by scattering off of heavy nuclei via a coherent 
scalar (“spin-independent”) interaction involving exchange of the lightest Higgs 
boson and the regions where capture occurs primarily by scattering via an axial 
(“spin-dependent”) interaction involving squark exchange off of hydrogen. We 
also show the regions of parameter space where the capture and annihiliation 
rates are large enough that the annihilation rate is half the capture rate and the 
neutrino flux is at “full signal.” 

In Section III we discuss the neutrino spectra ( dN/dE)fi from products 
of neutralino-neutralino annihilation in the Sun and Earth. We describe the 
hadronization and decays of the annihilition products and the interaction of the 
annihilation products and high-energy neutrinos with the Sun. In Section IV we 
discuss detection of high-energy neutrinos from the Sun (and Earth) and argue 
that for heavy WIMPs the neutrino signal from the Sim should be stronger than 
that from the Earth. We then point out that the most promising method of 
detection is via observation of upward-moving throughgoing muons induced by 
high-energy neutrinos in the rock below the detector and discuss the calculation 
of the event rate. 

In Section V we present our results, discuss which supersymmetric candidates 
for the primary component of the galactic halo are already ruled out by current 
neutrino-flux limits and which may be observable in the near future. Most of the 
models that axe inconsistent with current limits from IMB and Kamiokande 
on high-energy neutrino fluxes are those where the neutralino is a mixed gaug- 
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ino/Higgsino state and the mass of the lightest Higgs boson is near the current 
lower limits imposed by LEP. We find that if observational neutrino-flux limits 
are improved by a factor of ten, say, many more supersymmetric models will be- 
come detectable by these methods. The neutrino signal from neutralinos that are 
primarily gaugino is greater for models where the squark mass is smaller, while 
the neutrino rates from neutralinos that are Higgsinos or mixed gaugino/Higgsino 
states are relatively insensitive to the squark mass. In the last Section we dis- 
cuss our results, briefly discuss backgrounds and detection strategies, and make 
some concluding remarks. In Appendix A we display the cross section for elastic 
scattering of a neutralino off of nuclei, and Appendix B contains new results for 
cross sections for annihilation of neutralinos into mixed Higgs/gauge boson final 
states. 


II. RATE OF ANNIHILATION IN THE SUN 

The first step in calculating the rate for WIMP-induced neutrino events from 
the Sun is the determination of the rate at which neutralinos annihilate in the 
Sun. As mentioned previously, neutralinos accumulate in the Sun or Earth by 
capture from the galactic halo and are depleted by annihilation. If N is the 
number of neutralinos in the Sun then the differential equation governing the 
time evolution of N is 

N = C - C a N 2 , (4) 

where the dot denotes differentiation with respect to time. Here, C is the rate 
of accretion of neutralinos onto the Sun (or Earth). The determination of C is 
straightforward and will be discussed in detail below, and if the halo density of 
neutralinos remains constant in time, C is of course time-independent. 

The second term on the right-hand side is twice the annihilation rate in the 

Sun (or Earth), = C A N 2 /2, and accounts for depletion of neutralinos. The 

quantity C A depends on the cross section for neutralino-neutralino annihilation 

30 

and the distribution of neutralinos in the Sun (or Earth), 




(av) A V 2 

V? ' 


( 5 ) 
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where ( crv) A is the spin-averaged total annihilation cross section times relative 

velocity in the limit of zero relative velocity (since captured neutralinos move 

very slowly), and can be evaluated using the formulas in GKT and Appendix B, 

30 17 

and the quantities Vj are effective volumes for the Sun or Earth: 



where T is the temperature of the Sun or Earth, mpj is the Planck mass, and 
p is the core density of the Sun or Earth. In Ref. 30 it is found that Vj = 
6.5 x 10 28 (;'ml°) -3 / 2 cm 3 , where ml 0 is the neutralino mass in units of 10 GeV, 
for the Sun, and in Ref. 17 it is found that Vj - 2.0 x 10 25 (;'mi°) _3/2 cm 3 for the 
Earth. 

Solving Eq. (4) for N, we find that the annihilation rate at any given time is 

T A = j tanh 2 (t/r A ), (7) 

where t a = [CC A )~ X ^ 2 is the time scale for capture and annihilation to equili- 
brate. Therefore, if the the age of the Sun is much greater than the equilibration 
time scale (f 0 = 1.5 x 10 17 s » t a ) then the neutrino flux is at “full signal” 
(T a = C/2), but if Tj 4 > t 0 then the annihilation rate is smaller and the neu- 
trino signal is diluted accordingly. As we shall see, the capture rate in the Earth 
is generally < 10 -9 that in the Sun while the value of Vj in the Earth is only 
about 3 x 10 -4 that in the Sun, so the value of t a is always larger in the Earth 
than in the Sim; consequently, the fraction of full signal in the Earth can never 
be greater than that in the Sun. 

Although the calculation of the rate of accretion of WIMPs onto an astro- 
physical object is quite involved the basic idea is simple. Suppose a halo WIMP 
which has a velocity v far away from the object has a trajectory that passes 
through the object. At a point within the body where the escape velocity is Uesc 
the WIMP velocity will then be (w 2 , + u 2 *.) 1 / 2 . K the WIMP elastically scatters 
off of a nucleus of mass m, to a velocity less than u esc the WIMP will be captured. 
Kinematics tells us that the fractional energy loss (A E/E) of the WIMP in the 


9 


collision must lie in the range 


„ ^ A E . 4m cm, 
0<~ <— f— -J, 
Jb { m x + m i) 


( 8 ) 


and in the simplest case the cross section Og D for elastic scattering of the neu- 
tralino off of nucleus i is isotropic so the probability for a given energy loss is flat 
in this interval. [As will be discussed below, if the neutralino interacts coherently 
with the entire nucleus, at high momentum transfer there will be a form-factor 
suppression to the cross section so the probability for a given energy loss will 
no longer be flat in the interval given by Eq. (8).] The rate of capture of the 
WIMP by scattering off of nucleus i at this point in the Sim is then the rate of 
elastic scattering *SD n '(®«c + t &) 1/2 (where n' is the number density of nucleus 
i ) times the consitional probability that the WIMP is scattered to a velocity less 
than v esc : 



where x± = 4 m^m,/(m^ ± m,) 2 and 9 is the Heaviside step function. 

The essence of Gould’s resonant enhancement in the capture of WIMPs [and 
the kinematic suppression factor Si(m^) discussed below] is contained in Eq. (9): 
The conditional probability that a WIMP will be captured in a scattering event 
is greatest when is maximized which occurs when the neutralino mass closely 
matches the mass of the nucleus off of which it scatters. Furthermore, this 
resonance effect is much sharper in the Earth than in the Sun: The velocities 
of the WIMP have a Maxwell-Boltzmann distribution with velocity dispersion 
of v = 300 km s -1 and the escape velocity from the Earth ranges from 11.2 km 
s -1 (at the surface) to 14.8 km s -1 (at the center), so the probability is nonzero 
only for the very slow WIMPs on the Boltzmann tail or for WIMPs with masses 
that very nearly match m,. In a detailed analysis Gould finds that WIMPs 
in the “resonance range” 10-75 GeV have masses which are sufficiently close to 
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the mass of an element with a significant abundance in the Earth so that their 
capture is not kinematically suppressed. On the other hand, the escape velocity 
just at the surface of the Sun is 618 km s -1 (and Vesc is much greater at the 
center), so capture is not kinematically suppressed unless is quite small (i.e., 
the neutralino and nuclear masses axe very mismatched) and the resonance range 
for capture in the Sun is much larger than in the Earth. 

The neutralino scatters off of nuclei with spin (which for the purpose of 

capture in the Sun of Earth includes only the hydrogen in the Sun) via an axial 

or “spin-dependent” interaction characteristic of Majorana particles. In addition, 

the neutralino may scatter off of any nucleus via a scalar interaction in which 

the neutralino interacts coherently with the entire nucleus; for heavy neutralinos, 

the scalar cross section <xsc Is proportional to the fourth power of the nuclear 

5 31 

mass. For the elastic scattering cross section we use the results of Griest, 
which include both a spin-dependent and a scalar term due to the exchange of a 
squark and the Z boson, and of Barbieri, Frigeni, and Roulet, which includes a 
coherent scattering term due to the exchange of the lightest Higgs boson. We also 
include the effect of the exchange of H^, the heavier scalar Higgs boson (which 
increases the elastic scattering cross section only slightly). As recently pointed 
out by Gelmini, Gondolo, and Roulet, 21 the cross section for scalar interactions 
of neutralinos with nuclei is larger than that given in Refs. 5 and 32 when one 
takes into account the substantial strange-quark content in the nucleus as implied 
by the pion-nucleon sigma term. 3 For the convenience of the reader the complete 
formulas for the elastic scattering cross section are listed in Appendix A. 

Until now we have assumed that the elastic scattering cross section is isotropic 
and the conditional probability for a given energy loss in the range given by 
Eq. (8) is uniform; however, if the neutralino interacts coherently with the nucleus 
and the momentum transfer q is not small compared to the inverse of the nuclear 
radius R this assumption is not necessarily true as the neutralino does not see 
the entire nucleus and the cross section for scattering of neutralinos off of nuclei is 
form-factor suppressed (like that for electromagnetic elastic scattering of electrons 
from nuclei). In terms of the energy loss AE the form factor suppression may be 
written as 34 
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[F(9 2 )| 2 = exp(-A E/Eo) 


( 10 ) 


where Eq = 3/(2 rriiR 2 ). 

Now let us consider the relevance of a form-factor suppression on the capture 
of heavy neutralinos in the Sun and Earth. First of all, for a WIMP with a kinetic 
energy = m^v^/2 in the halo to be captured it must have an energy loss in 
the range 

£oo — A E < X+(E<x> 4“ ^esc)i (11) 

where E e sc = m^v“ sc /2 is the WIMP escape energy at the point of collision in 
the Sun. The lower limit comes from the condition that the WIMP scatter from 
a velocity ( v ^ to a velocity less than Vesc, and the upper limit is the 

kinematic limit. This implies that in order to be captured the WIMP energy in 
the halo must be E 0 0 < x~E e sc, which in turn implies that the largest energy loss 
involved in capture of WIMPs from the halo is AE miX = X-Eesc- The value of Eq 
for iron, the heaviest element important for capture in both the Sun and Earth, 
is 8 x 10 -5 GeV. Because of the factor of (m^ — m,) 2 in the denominator of X- 
the energy loss is largest for the lightest WIMP we consider, one with a mass of 
80 GeV. For capture in Earth, the largest energy loss occurs at the center of the 
Earth and is roughly 2 x 10 -6 GeV, so the form-factor suppression is negligible 
for capture of heavy WIMPs in the Earth. On the other hand, the maximum 
energy loss for capture off of iron in the Sun is 8.1 x 10 -3 , which implies that a 
proper calculation of capture in the Sun must include the effects of form-factor 
suppression of the coherent scalar interaction. 

The full capture rate calculation assumes the astrophysical object moves 
through a homogeneous Maxwell-Boltzmann distribution of WIMPs and requires 
information about the elemental composition of the object and the distribution 
of elements in the object. One must integrate over the trajectories of the WIMP 
through the Sun and over the velocity distribution of the WIMPs. The fined 
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result for the capture rate, adapted from Gould, is 


C = 


. Poa 
m^v 300 


[ <7 SD° ) + F «( m *) cr sc 0) 


( 12 ) 
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where c = 5.8 x 10 24 s -1 for the Sun and c = 5.7 x 10 15 s -1 for the Earth, is 
the mass density of neutralinos in the galactic halo in units of 0.4 GeV cm -3 , 
is the neutralino mass in units of GeV, and U 300 is the velocity dispersion of the 
neutralinos in the galactic halo in units of 300 km s _1 . The sum is over all species 
of nuclei in the astrohysical object (here the Earth or Sun), mi is the mass of 
the zth nuclear species in GeV, /i is the mass fraction of element z, cTgp 0 ^ is the 
cross section for elastic scattering off of nucleus i via an axial interaction (given 
in Appendix A) in units of 10 -40 cm 2 , and <Xg^°* is the cross section for elastic 
scattering of the neutralino off of nucleus z via a coherent scalar interaction (given 
in Appendix A) in units of 10 -40 cm 2 . The quantities <f>i describe the velocity 
distribution of element z in the Sun or Earth and are given in the Appendix of 
Ref. 21 as are the quantities /,. 

The quantity 5,(m^) is the kinematic suppression factor for capture of a 
WIMP of mass off of nucleus z. We use an approximation that interpolates 
between the two limiting cases of the suppression factor given by Gould: 




A b 

1 + A b 



(13) 


where 


and b = 1.5. We obtain this expression from the RS expression which approxi- 
mates Gould’s kinematic factor to 5% for scattering off of protons by noting that 
the neutralino and nuclear masses enter into Gould’s kinematic formula only in 
the combination m^m,/(m^ — m,) 2 . The quantity t>esc is the escape velocity at 
the surface of the Sun or the Earth (618 km s -1 for the Sun and 11.2 km s 1 
for the Earth). To check our approximation for Si , we calculate the capture rate 
in the Earth for neutralinos with masses between 10 GeV and 80 GeV using an 
elastic scattering cross section due only to neutral-Higgs and Z exchange and 
note that our results reproduce the resonance structure found in Ref. 20. Note 
also that Eq. (12) reduces to the simple expression for capture by the Earth when 
the mass of the neutralino is far from the resonance range (Ref. 18). 


me — miY \ v J 


(14) 


13 



Although heavy neutralinos are outside the resonance range 10 GeV < < 

75 GeV for capture by the Earth so that only the form of Si in the limit A < 1 
is important for capture by the Earth, is much larger for the Sun and is 
typically 2-3 times larger for the Sun, so as mentioned above, the resonance range 
for capture in the Sun is much wider than in the Earth. For example, in the Sun 
A = 1 when the neutralino mass is roughly seven times that of the nucleus off of 
which it scatters. 

The form-factor suppression F,(m^ ) of the capture of a WIMP of mass 
from nuclei i is obtained by comparing the results of integrating Gould’s differ- 
ential capture rate [Eq. (A10) in Ref. 17] over the mass of the Sun including 
form-factor suppression and comparing it with the integral of the analogous ex- 
pression [Eq. (2.24) in Ref. 17] in which full coherence is assumed. In doing so, 
the density of the Sun as a function of radius r was taken to be 

p(r) = po exp(-7.7r/jR 0 )(l - r/R 0 ) 16 , (15) 

where po is the density at the center of the Sun, and R@ is the solar radius. This 
form approximates the solar density in Ref. 35 and yields the correct gravitational 
potential at the center of the Sun (5.1 times as large as the potential at the surface 
of the Sun) and the average gravitational potential for heavy nuclei (3.4 times 
that at the surface of the Sun). The resulting Fi are plotted in Fig. 2. From 
Fig. 2 we see that the form-factor suppression for capture from scattering off of 
hydrogen and helium is negligible, capture from scattering off of elements with 
atomic masses 12-32 is moderately suppressed, while capture from scattering off 
of iron is suppressed by several orders of magnitude for WIMPs in the several 
hundred GeV range. If there were no form-factor suppression, owing to the 
factor of m* in the scalar cross section one would expect scattering from iron 
nuclei to dominate the capture of WIMPs in the Sun; however, because of the 
form-factor suppression, capture of heavy WIMPs in the Sun occurs primarily 
by scattering off of oxygen. Even so, capture from scattering off of iron nuclei 
is still significant. When considering the complete capture rate due to scalar 
interaction of WIMPs off of nuclei in the Sun, one finds that the form-factor 
suppression of the scalar elastic scattering cross section decreases the capture 
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rate by a factor of about 0.3 for WIMPs of mass 80 GeV and about 0.07 for 
TeV-mass WIMPs. 

The relative importance of the capture rates due to spin- dependent scattering 
opposed to coherent scattering due to squark and Higgs exchange depends on the 
supersymmetric model. Coherent scattering vanishes as the neutralino becomes 
a pure B-ino or Higgsino as does spin-dependent scattering due to Z exchange. 
To study the effect of Higgs-exchange scattering on the capture rate we set the 
squaxk mass to infinity. Doing so we find that the capture rate due to Higgs 
exchange is generally more important than that due to Z exchange when the 
neutralino is heavier than the W. In Fig. 3 we show contour plots in the M-/x 
plane of the rate of capture of neutralinos in the Sim for (a) tan /3 = 2, m H o = 20 
GeV, and n > 0; (b) tan/3 = 2, m H o = 20 GeV, and fi < 0; (c) tan/3 = 2, 
m H o = 35 GeV and n > 0] and (d) tan/3 = 25, m H o = 35, and n > 0 assuming 
the squaxk mass is infinite. As expected, when squark exchange is negligible 
mixed-state neutralinos Eire captured far more readily than pure B-inos or pure 
Higgsinos. For fixed masses the capture rate decreases with increasing purity. For 
heavy neutralinos of fixed gaugino/higgsino composition that are heavy enough 
to be outside the Sun’s resonance range, the capture rate is roughly proportional 
to m* -5 / 2 ; one factor of m* -1 is due to the number density in the galactic, one 
factor of m^ -1 is due to the kinematic suppression [c/., Eqs. (13) and (14)], while 
the other factor of roughly m^ -1 / 2 comes from form-factor suppression. [The 
cross section for scattering due to exchange of the lightest Higgs boson does not 
decrease as the neutralino mass is increased fax past the mass of the nucleus off 
of which it scatters; see Eq. (2.5) in Ref. 20.] Incidentally, as the neutralino 
mass is increased past a TeV, the form-factor suppression ceases to decrease 
with increasing WIMP mass; the reason is that if the nuclear mass is negligible 
compared to the WIMP mass, the momentum transfer does not depend on the 
WIMP mass. 

From Fig 3 we also find that if tan /3 is held fixed, the capture rate generally 
decreases with increasing m go due to the propagator suppression, and if we 
hold m H o fixed, the capture rate generally increases with increasing tan 0; this is 
simply because the Higgs couplings contain terms inversely proportional to cos/3. 
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To see the 'effect of the squark mass on the capture rate we show in Fig. 4 the 
rate of capture of neutralinos in the Sun when we take the squaxk mass to be 20 
GeV heavier than the neutralino mass. Doing so, we find that the capture rate for 
Higgsinos and mixed-state neutralinos is similar to that when the squark mass is 
infinite; this implies that capture of Higgsinos and mixed-state neutralinos occurs 
primarily by Higgs-exchange scattering and the capture rate is insensitive to the 
squark mass. On the other hand, for models where the neutralino is mostly B-ino 
and the squaxk is taken to be 20 GeV heavier than the neutralino, capture occurs 
primarily by spin-dependent scattering of the neutralino off of the hydrogen in 
the Sun. This is illustrated in Fig. 5 where we show contours of the fraction of 
the capture rate that occurs due to spin-dependent scattering. Scattering that 
occurs via spin-dependent exchange of the squark depends only very weakly on 
tan /? and does not depend on m y o at all; therefore, if the squaxk mass is small 
enough so that capture of the neutralino occurs primarily by squark- exchange 
scattering, the capture rate depends primarily on the squaxk mass. We should 
also mention that in computing the spin-dependent cross section we used the (still 
controversial) EMC results for the spin content of the proton. As discussed 
in Appendix A, if instead we used the naive flavor-SU(3) quark model for the 
proton the spin-dependent cross section due to squark exchange would be rougly 
3 times larger. 

Now that we have results for the capture rate we can see where the annihila- 
tion rate is at full signal Tyt = C/2 and where the time scale for equilibration of 
the number of WIMPs N is so large that < C. In Fig. 6 we show the regions 
of parameter space where energetic neutrinos are not at full signal because neu- 
tralinos have not had sufficient time to collect in the Sun. In the dark shaded 
regions the signal is less than 10% of the full signal (f©/rA < 0.33) and in the 
light shaded region the signal is l6ss than 90% of the full signal (to/r A < 1.82); 
elsewhere, capture and annihilation of neutralinos occurs rapidly enough so that 
the neutrino rates are at full signal (to/r^ > 1.82). In Fig. 6(a) tan/? = 2, 
mjyo = 20, the squark mass is taken to be infinite, and fi > 0; Fig. 6(b) is similar 
but n < 0 is shown; and Fig. 6(c) is similar to Fig. 6(a) but the squaxk mass is 
taken to be 20 GeV heavier than the neutralino mass. Note that in most models 
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where the neutralino is lighter than a TeV the neutrino flux is at full signal. 
Later we will find that in regions of parameter space where the neutrino flux is 
laxge enough to be near current observational limits, the flux is at full signal. 
We will also see that t\ generally stays small enough so that the rates remain at 
full signed even for most models with a neutrino flux severed orders of magnitude 
weaker than the current observational limits. 

III. NEUTRINO SPECTRA FROM NEUTRALINO ANNIHILATION 

Given the annihilation rate Ta = C tanh 2 (*©/fA)/2, the differential flux of 
neutrino type i (e.g., u e , etc.) produced by the annihilation of neutralinos 

in the Sun or Earth at a distance R from the source is 



where the sum is over all annihilation channels. 

The quantities Bp are the branching ratios for annihilation into final state F. 
Since the neutralinos are moving nonrelativistically in the Sun or Earth, Bp may 
be determined by the relative magnitude of the cross sections for annihilation into 
channel F at zero velocity given in Ref. 6 [Eqs. (A10), (B7), (Cll), and (D6)] 
and in Appendix B. The final states F into which neutralinos may annihilate at 
zero relative velocity are ff where / is a quark or charged lepton, W + W , Z°Z°, 
Z°H%, Z°H%, W+H~, and W~H + 37 The cross sections for annihilation 
into other combinations of gauge and Higgs bosons vanish as the relative velocity 
approaches zero. The calculation of the cross sections for annihilation into the 
mixed gauge- and Higgs-boson final states W+H , and W H+ at 

zero relative velocity are new and the results are presented in Appendix B. As 
noted by Olive and SrednickiJ annihilation into the mixed gauge/Higgs boson 
final states is generally small for pure iMnos and Higgsinos but may be important 
for mixed-state neutralinos. For models where the neutralino is a pure U-ino and 
the squark masses are much larger than all other masses involved, annihilation 
into the mixed gauge/Higgs boson states may be comparable to annihilation into 
Higgs-boson states; in this case, neutralinos annihilate predominantly into these 
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states, but the total rate for annihilation is very small and the neutralinos are 
generally very weakly interacting. 

The ( dN/dE)fi are the differential energy fluxes of neutrino type i at the 
surface of the Sun (or Earth) that result from the injection of particles in final 
state F at the center of the Sun (or Earth). These fluxes are functions of the 
neutrino energy E and the energy E, of the injected particles. Calculation of 
the fluxes requires information about the cascade following the decay of the 
annihilation products, the hadronization of heavy quarks in the cascade, and 
the interactions of particles in the cascade with the medium at the core of the 
Sun or Earth. Since the ( dN/dE)n are the neutrino fluxes at the surface of 
the astrophysical object while neutralino annihilation occurs at the center and 
the Sun is not transparent to neutrinos with energies in the 100-GeV range, 
absorption and energy loss of neutrinos by the solar medium must also be included 
in the calculation. Since the density and thickness of the Earth Eire different 
from those in the Sun, the ( dN / dE)f{ from particles injected in the Earth will 
be different than those from the Sun. 

A detailed calculation of the neutrino spectrum from injected quarks and lep- 
tons was performed by the authors of Ref. 19 using the Lund Monte Carlo. Their 
calculation includes hadronization of quarks and interactions of the fermions and 
neutrinos with the solar medium. Electrons, muons, and light (u, d, s) quarks are 
stopped in the Sun before they decay and therefore do not produce high-energy 
neutrinos. The top quark is expected to hadronize and then decay far before it 
can lose a substantial amount of energy, and the r will also decay immediately. 
Bottom and charm quarks hadronize and due to the high density of the core 
of the Sim, the heavy hadron may subsequently lose a significant fraction of its 
energy before decaying. RS estimate that if Eq is the initial heavy hadron energy 
in the Sun, the mean energy of the hadron when it decays will be 

00 

(E) = E c e Ec ! Eo J ^dx, (17) 

Ec/Eo 

and they estimate that E c ~ 250 GeV for charmed hadrons, and E c « 470 GeV 
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for bottom hadrons. Evaluating the integral, one finds that (E) ~ Eo(l — Eo/E c ) 
for Eo < E c and (E) ~ E c [\n(E 0 /E c ) - 7 e] for E 0 > E c , where 7 e = 0.577... is 
Euler’s constant, so the mean energy of the decaying hadron never grows much 
larger than E c . 

At high energies the Sun is no longer transparent to neutrinos and interactions 
of neutrinos with the solar medium may significantly alter the energy spectrum. 
For r’s injected at energies above several hundred GeV, the flux of muon neu- 
trinos may be significantly enhanced by the decay of additional r’s produced by 
.charged-current interactions of tau neutrinos with the solar medium. Electron 
and muon neutrinos are absorbed by charged-current interactions: The probabil- 
ity that a neutrino of initial energy E, will escape from the Sun is exp(— E%J E Ab s)i 
where E Abs = 198 GeV for neutrinos and E Abs = 296 GeV for anti-neutrinos. 
Furthermore, at high energies neutral-current interactions degrade the neutrino 
energy. 

Since the density of the core of the Earth is about 1/12 that of the core of the 
Sun, muons and light quarks are still stopped before they decay, while stopping 
of heavy hadrons may be ignored until several TeV. Moreover, the optical depth 
of the Earth is much smaller than that of the Sun, so interactions of neutrinos 
with the Earth may be ignored for neutrino energies less than several TeV. As 
a result, Ritz and Seckel’s non-interacting results may be used for the neutrino 
spectra from the Earth. 

The results presented by RS are for neutrino spectra from fermions injected 
into the core of the Sun at 60 GeV and at 1000 GeV; however, we need to 
obtain information about the spectra for fermions injected at any energy up to 
a TeV. For re as ons to be discussed below, we will eventually focus on detection 
of neutrinos via neutrino-induced upward-moving muons. Since the cross section 
for a neutrino to produce a muon in the rock below the detector is proportional 
to the neutrino energy and the range of the muon is roughly proportional to 
the energy, the probability for a neutrino to produce a throughgoing muon is 
proportional to the energy squared. Therefore, to obtain event rates we need 
only the second moments (N z 2 ) ^m^ 2 , where 
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E 2 dE. 


(18) 


The functional forms of the spectra are not required. 

For fermions injected into the core of the Earth, interactions are negligible and 
the moments of the neutrino spectra are easily obtained from Ritz and Seckel’s 
non-interacting results. In this case, 19 

{Nz 2 ) = \{N){y 2 ){{z})-^j, (19) 

where E t is the fermion injection energy, (TV), (y 2 ), are the rest-frame yield 
and second moments listed in Table 2 of RS, is the second moment of the 

fragmentation function listed in Table 3 of RS, and my is the mass of the injected 
fermion. 

For fermions injected into the core of the Sun, the calculation is much more 
difficult since one must take interactions into account. RS outline a procedure for 
analytically estimating the the effect of interactions which reproduces the Monte 
Carlo results reliably for injection energies < 200 GeV. An effort to modify 
and apply the corrections to describe interactions at higher energies resulted in 
moments of the neutrino spectra that reproduced those obtained from the Monte 
Carlo only to within ~ 50%; however, in doing so one finds that for injected b and 
c quarks, the most important effect is the stopping of heavy hadrons. Therefore, 
for the scaled second moment of the neutrino spectra for b and c quarks we 
assumed that 

2 

, ( 20 ) 

where xo = E c /Ei, and fitted a and E c to match the interacting results of RS 
at 60 GeV and 1000 GeV. (Actually, since RS did not present interacting results 
at 60 GeV for anti-neutrinos or, in the case of the b quark, for neutrinos, we 
obtained these numbers using the corrections for interactions described in their 
paper.) We found that for neutrinos from c quarks, a = 0.056 and E c = 155 GeV; 
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for anti-neutrinos from c quarks, a = 0.052 and E c = 275 GeV; for neutrinos from 
b quarks, a = 0.086 and E c = 185 GeV; and for anti-neutrinos from b quarks, 
a = 0.082 and E c = 275 GeV. 

Since r leptons are not stopped and do not hadronize, absorption of muon 
neutrinos is the most important interaction effect for the spectra from r leptons; 
production of muon neutrinos from interactions of r neutrinos is also significant 
at high energies, but these neutrinos axe predominantly low energy and do not 
contribute significantly to the second moment. Thus, we take the second moment 
of the neutrino spectrum from injected r leptons to be 

(Nz 2 ) = ae~ E ' /EKb ‘, (21) 

and fit a and E\\, s to reproduce the RS results at 60 GeV and 1000 GeV. For 
neutrinos, a = 0.0204 and £ A bs = 476 GeV, and for anti-neutrinos, a = 0.0223 
and jE/Abs = 599 GeV. 

Our estimates for the spectra from the top quark are far more uncertain. RS 
used a top-quark mass of 40 GeV, and here we have assumed that it is 120 GeV. 
Since even 40 GeV is so much heavier than all other lighter particle masses, we 
assumed that the scaled rest-frame neutrino spectra would be the same for a top 
quark of 120 GeV as it would for a top quark of 40 GeV. We then estimated the 
effect of interactions for a top quark injected into the solar core at 120 GeV and 
assumed that the RS interacting results at 1000 GeV would also be valid for a 120 
GeV top quark. At injection energies just above threshold the moments of the 
neutrino distribution have a strong dependence on the fragmentation function, 
and at higher energies, absorption of neutrinos determines the behavior of the 
spectral moments. Therefore, neither of the expressions in Eq. (20) or Eq. (21) 
really describe the injection-energy dependence of (Nz 2 ). Nevertheless, the effect 
of interactions, which we can reliably estimate at low energies, is better described 
by Eq. (20) than by Eq. (21), so we use the form of Eq. (20) with a = 0.18 and 
E c = 110 GeV for neutrinos, and a = 0.14 and E c = 380 GeV for anti-neutrinos. 

Although these estimates of (Nz 2 ) are somewhat ad hoc and admittedly 
crude for arbitrary injection energies between 60 GeV and 1000 GeV , they should 
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be relatively accurate for neutrino spectra from annihilation of neutralinos not 
much heavier than the W or Z; at higher energies, our approximations are far 
from pinpoint accuracy, but they should still be good enough to indicate the 
effect of interactions of the decay products and neutrinos with the solar medium. 

Since Higgs and vector bosons decay into pairs of quarks and leptons imme- 
diately, it is easy to obtain (N z 2 ) for injected bosons from our previous results 

.38 

for the neutrino spectra from injected fermions. Suppose boson B undergoes 
2-body decays into fermions /, and N j is the number of fermions of type / pro- 
duced on average per B decay (i.e., ^ y = 2) and can be obtained from the 
branching ratios for decay of B into the various final states and the contents of 
those channels. If E, is the injected boson energy, then the energy of the fermion 
in the rest frame of the B is mg/ 2, where mg is the B mass, and in the moving 
frame it is Ef = E,( 1 -f /?cos#)/2 where 0 is the velocity of B in units of the 
speed of light and 6 is the angle between the direction of motion of the decay 
product and the direction of motion of B. For Higgs bosons and unpolarized 
vector bosons (which are produced by the annihilation of neutralinos provided 
the interactions of the neutralinos are CP-conserving, which is assumed through- 
out here) the decay is isotropic which means that the laboratory-frame energies 
of the fermions from the decay of B are evenly distributed from E t (l — ft)/ 2 to 
Ei(l + 0)/2. Therefore, 


B E,( l+0)/2 

( n * 2 )b. = Zw / < 22 > 

f '«< i-m 

where i^Nz 2 ') ^ (E) are the second moments of the neutrino spectra presented 
above as a function of the injected fermion energy E. 

The three neutral Higgs bosons of the minimal extension of the supersymmet- 
ric standard model decay into fermion- antifermion pairs. The branching ratios 
for the decays of ^d H® , from which the Nf are obtained are given in the 
Appendix of Ref. 21 and are proportional to the fermion mass squared (so the 
Higgs bosons do not decay directly into energetic neutrinos), and the branching 
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ratios for the decay of H\ may be obtained from those for if® decay by switching 
cos a and sin a. 

If the neutralinos annihilate into r leptons or 6, c, or t quarks, and an energetic 
neutrino is produced in the decay of these fermions, then the typical neutrino 
energy is 1/3 the mass of the neutralino. If the neutralinos annihilate into Higgs 
bosons there is another step in the decay chain before energetic neutrinos axe 
produced, so their energies would typically be 1/6 the mass of the neutralino. 
This is partially compensated by the fact that each Higgs boson produces two 
fermions, but since the detection rates are proportional to the energy squared, the 
net effect is that if the neutralinos annihilate into Higgs bosons the detection rate 
is roughly half the rate if they annihilated into fermions (assuming, of course, 
that the branching ratio for the various fermions from Higgs-boson decays is 
nearly the same as the branching ratios for the various fermions from neutralino 
annihilation if only fermion final states are considered). Although H® must be 
lighter than mg cos 2/3, and most certainly decays only into quarks and leptons, 
the other Higgs bosons may be much heavier and may include other exotic decay 
channels as well which may also produce energetic neutrinos which would most 
likely have a much softer spectrum. If this is the case, then by assuming that they 
decay only into quarks and leptons we are overestimating the neutrino yields. 

It turns out that the most favorable annihilation channel for observing high- 
energy neutrinos is the gauge-boson fined state. The reason is that W and Z 
bosons decay directly into neutrinos with appreciable branching ratios. Com- 
pared with the event rate from these “semi-prompt” neutrinos, the event rate for 
neutrinos which come from the quark and charged-lepton decay products of the 
gauge bosons is negligible. A W decays to a muon and a muon neutrino about 
11% of the time, 39 so neglecting interactions, { Nz 2 ) is roughly 0.025 for slow _ 
W’s and 0.033 for relativistic W’s. This is larger than all the values expected 
from fermion- antifermion pairs (see Table 1 of RS), although final states come 
close. Furthermore, at higher energies, no energy is lost from hadronization or 
stopping of the vector bosons. (At higher energies, the value of (IV z 2 ) for gauge- 
boson final states becomes smaller than that from r ± final states; this is because 
the energies of neutrinos from gauge-boson decays are generally larger than those 
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from t decays so absorption of neutrinos in the Sun from gauge-boson decays is 
stronger than absorption of neutrinos from r decays. Even so, if the neutralino 
annihilates to r* pairs, it will also have a significant and usually larger anni- 
hilation branch to 66, cc, and if kinematically accessible, tt pairs, so the total 
neutrino yield from gauge-boson final states will be greater than the total yield 
from fermion-antifermion states.) The branching ratio for Z° — ► vv is slightly 
smaller than the branching ratio for W —* but two neutrinos are produced 
so ( Nz 2 ) is a little larger. 

For W bosons injected in the core of the Earth with velocity we can ignore 
interactions of the neutrinos with the Earth, and 

( Nz2 )wi = IV— A«iv( 3 + /3 2 )/l2, (23) 


where i is a neutrino or anti- neutrino; (Nz 2 ) Zi may be obtained by multiplying 
by two and replacing by Tz^pp- To account for interactions of the 

neutrinos with the solar medium for vector bosons injected into the core of the 
Sun we use the estimate of RS that a neutrino injected with an an energy E 
leaves the Sun with energy 



E 

1 + Et{ ' 


(24) 


where t„ = 1.01 x 10 3 GeV 1 and Tp = 3.8 x 10 4 GeV 1 , and probability 


Pf = 



(25) 


where a„ = 5.1 and otp = 9.0 for anti-neutrinos. Doing so we find that 

,£=£.( l-0)/2 


/ pr 2\ _ IV— 1 2 ~t~ 2-E l T|( 1 4~ Qt) + + 0)') 

\ Iwi a,r 3 (a 3 — 1)(1 + Et,)®'* 1 


(26) 


E=E,( l+0)/2 


for W's injected into the core of the sun with energy E t . 

In Fig. 7 we show the second moments m% 2 (Nz 2 } of the neutrino yield from 
the Sun for the cc, 66, tt, r ± , W ± , and (using tan/3 = 2 and m#o = 20) 

final states as a function of the neutralino mass. The neutrino yields from Z° 
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pairs (not shown) is similar to, but slightly smaller, than the yields from W ± pairs 
and the yield from the (when it is kinematically accessible) final state is 

similar to that from the final state. We remind the reader that although 

the yield from r ± pairs surpasses that from gauge-boson pairs for neutralinos 
heavier than about 200 GeV, if the neutralino annihilates to lepton pairs, then it 
also has a significant annihilation branch into quark-antiquark pairs and the yield 
from gauge-boson pairs is still larger than the total yield from fermion- antifermion 
final states. 

IV. RATES FOR DETECTION IN UNDERGROUND DETECTORS 

Generally, neutrinos axe detected either by contained events where the neu- 
trino undergoes a charged- current interaction and produces a lepton in the detec- 
tor or by upward-moving throughgoing muons in which a muon neutrino under- 
goes a charged- current interaction in the rock below the detector and produces 
a muon which then passes through the detector. Since the cross section for a 
charged- current interaction is proportional to the neutrino energy and the effec- 
tive range of a muon is proportional to the muon energy, the rate for contained 
events is roughly proportional to the neutrino energy and the rate for neutrino- 
induced throughgoing muons is proportional to the square of the neutrino energy. 
Therefore, at sufficiently high energies the rate for throughgoing muons should be 
greater than that for contained events. In Ref. 21 the regions of parameter space 
ruled out by searches for contained events from Frejus 40 very nearly matches 
those regions ruled out by searches for throughgoing muons from IMB for 
neutralinos less massive than the W. (In addition, NUSEX reports that limits 
on muons produced by neutrino interactions in the rock below the detector that 
stop inside the fiducial volume of the detector are in agreement with those from 
contained or throughgoing events from IMB, Frejus, and Kamiokande.) Further- 
more, Ref. 20 indicates that the rate for detecting high-energy neutrinos from the 
Sun via throughgoing muons per 100 m 2 becomes larger than that for contained 
events per kiloton for neutralinos heavier than roughly 60 GeV while the rate 
for observing throughgoing muons is greater than that for contained events from 
neutrinos from the Earth for neutralinos heavier than roughly 20 GeV . Therefore, 
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since neutralinos heavier than the W are considered here, we will concentrate on 
detection of neutrinos via throughgoing muon events. 

After taking the cross section for muon production in the rock and the effec- 
tive range of the muons into account but ignoring detector thresholds (which are 

near 2 GeV — far lower than the average neutrino energies considered here), the 

.19 

rate (per unit detector area) for neutrino-induced throughgoing muon events is 

r*, detect = 1-27 X 10 _29 Cm S 2 ^ a * b ' Bf ( N Fi m ~ 2 y r_1 > ( 27 ) 

i F 

for neutrinos from the Sun; the same expression multiplied by 5.6 x 10 8 (the 
square of the ratio of the Earth-Sun distance to the Earth’s radius) gives the 
rate for neutrino events from the Earth. Here, C is the capture rate in units of 
s -1 , the sum on i is over muon neutrinos and anti-neutrinos, the a, are neutrino 
scattering coefficients, a v = 6.8 and ap = 3.1, the 6, are muon range coefficients, 
b v = 0.51 and bp = 0.67, and (iVz 2 )^ is the second moment of the spectrum of 
neutrino type i from final state F scaled by the neutralino mass squared. 

Given the expressions for (N z 2 ) for neutrino spectra from the Sun and the 
Earth it turns out that for heavy neutralinos with masses not much greater than 
a TeV, the neutrino signal from the Sun should be larger than that from the 
Earth. To see this, first note that the difference in the prefactors c for the Sim and 
the Earth in the capture- rate equation Eq. (12) is roughly compensated by the 
geometric factor 5.6 x 10 8 accounting for the difference in the distances between 
us and the Sun or the center of the Earth (5.8 x 10 24 s -1 for the Sim opposed 
to 3.2 x 10 24 s -1 for the Earth). Therefore, for dark-matter candidates with 
masses in the Earth’s resonance range 10 GeV < m* £ 75 GeV the kinematic 
suppression factor Sj is nearly unity and since the fraction of the Earth’s mass 
due to heavy elements is higher than that in the Sun, the neutrino flux from 
the Earth may well be comparable to or greater than that from the Sun (if the 
WIMP in question has spin- independent interactions). 

In contrast, the heavy neutralinos considered here have masses outside the 
Earth’s resonance range, so capture by the Earth is strongly suppressed due to 
the factor of (v esc /v) 2 « 1.4x 10 -3 . Even if the neutralino mass is large enough 
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to fall outside the Sun’s resonance range, capture of WIMPs by the Earth is 
still suppressed relative to capture in the Sun because v^c is so much smaller in 
the Earth than in the Sun. Although form-factor suppression does not occur for 
capture of heavy WIMPs in the Earth, the form-factor suppression of the capture 
rate in the Sun never falls far below 10 — 1 , whereas the kinematic suppression of 
capture of WIMPs in the Earth is of order 10 -4 that in the Sim. In addition, if 
the capture and annihilation rates for the neutralino in question are small then 
the neutrino signal from the Earth may be further weakened relative to that 
from the Sim as the time for the number of neutralinos to reach equilibrium 
in the Earth is generally smaller than that in the Sun. When considering heavy 
WIMPs the calculation of the capture rate in the Earth is also far more uncertain 
than that for capture in the Sun. The reason is that only heavy WIMPs that 
are moving very slowly may be captured in the Earth and an isotropic Maxwell- 
Boltzmann distribution does not necessarily give a good approximation to the 
phase-space density of such WIMPs. At this point we should also remind the 
reader that if the WIMP in question has only an axial interaction with nuclei 
(such as a B-ino in models with a relatively light squark) it may be captured in 
the Sun by scattering off of hydrogen, but it will not be captured in the Earth. 

Although the rate of capture of WIMPs by the Earth remains small relative 
to that by the Sun at higher WIMP masses, above some large mass the neutrino 
signal from the Earth might become comparable to or larger than that from the 
Sun because of interactions of decay products and neutrinos with the Sun. We 
can estimate this mass scale by taking the following simplified model: Assume 
capture occurs only by scattering off of iron and neglect form-factor suppression of 
capture in the Sun; in the Earth this provides a good estimate of the capture rate, 
and if anything, this should underestimate the capture rate in the Sun. Doing 
so, the capture rates differ only in the prefactors c, the factors /,, <f>i, and the 
factor of vl sc in 5,. To include the effect of form-factor suppression of capture 
in the Sun we multiply the capture rate by 0.07, the value of the suppression 
factor for a WIMP of mass 1 TeV. (At smaller WIMP masses the suppression is 
not as severe, while the suppression does not become significantly stronger for 
WIMP masses greater than a TeV.) Doing so we find that the capture rate of 
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very heavy neutralinos in the Sun is roughly 20 times that in the Earth when 
scaled by the difference in the geometric factor. It turns out that the values of 
(W 2 ) for particles injected into the core of the Sun at 1 TeV are of order 1/10 
of those for particles injected into the core of the Earth at 1 TeV; therefore, the 
neutrino signal from heavy neutralino annihilation in the Sun should remain much 
larger than that from the Earth for neutralino masses below 1 TeV and become 
comparable at a WIMP mass near 1 TeV. For WIMPs heavier than a TeV, the 
neutrino signal from WIMP annihilation in the Earth should become stronger 
them that from annihilation in the Sun because of absorption of neutrinos in the 
Sun. Here we have ignored the fact that the number of very heavy neutralinos 
in the Sun or Earth may not have reached equilibrium; correcting for this would 
only increase the mass scale at which the neutrino signal from the Earth might 
become comparable to that from the Sun. 

So, the neutrino signed from the Sun should be much stronger than that from 
the Earth for neutralinos just heavier than the W, and the strength of the signal 
from the Sun relative to that from the Earth should decrease as the neutralino 
mass is increased until a WIMP mass of order a TeV when the signal from the 
Earth becomes comparable to that from the Sun. Since the signal from the 
Earth should be small compared to that from the Sun in the range 80-1000 TeV, 
in the following we will focus our attention on the neutrino signed from WIMP 
annihilations in the Sun only. We should also point out that these results imply 
that observation of a neutrino signal from the Sun and the absence of one from 
the Earth would be a signature of particle dark matter in the mass range 80-1000 
GeV. 


V. RESULTS 

Since the MSSM has many undetermined parameters we will show results in 
the M-n plane for several values of tan 0 and m # o allowed by null results from 
searches for neutral Higgs bosons at LEP . Again, we will first take the squark 
masses to be infinite; this minimizes the capture rate and emphasizes gauge- and 
Higgs-boson final states. Then we will consider squark masses 20 GeV higher than 
the neutralino mass; this will emphasize capture by spin-dependent scattering and 
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fermion final states for neutralinos where such effects are important. 

When the neutralino is mostly Higgsino, it annihilates primarily into gauge 
bosons, and the effects of the squark, Higgs-boson, and top-quark masses axe 
relatively unimportant. When the neutralino is mostly B- ino, it annihilates 
primarily into fermions (provided the squark mass is not too large), and when 
the top-quark channel is open, it annihilates predominantly into the top quark. 
Mixed-state neutralinos generally annihilate into gauge bosons, fermions, and 
Higgs bosons as well with comparable magnitudes. 

In Fig. 8 we plot contours of the fraction of the neutrino signal that comes 
from gauge bosons. When the squark mass is taken to be infinite [Fig. 8(a)], the 
neutralino does not annihilate into fermions and since gauge bosons yield a much 
harder spectrum of neutrinos than Higgs bosons, virtually all of the neutrino 
signal from heavy neutralinos comes from gauge-boson final states. When the 
squark mass is 20 GeV heavier than the neutralino mass [Fig. 8(b)], fermions 
are the dominant annihilation products from H-inos, so the neutrino signal is 
not always dominated by neutrinos from gauge bosons. Still, neutrinos from 
gauge-boson final states dominate the signal for Higgsinos and contribute a signal 
comparable to that from fermions in many regions of parameter space with mixed- 
state neutralinos and H-inos. 

The IMB collaboration has found an upper limit on the flux of upward- 
moving muons induced by neutrinos from the Sim with energy larger than 2 GeV 
of 2.65 x 10 -2 m -2 yr -1 , 28 (and similar, though slightly weaker limits have 
been fo un d by Kamiokande II 29 ). Therefore, supersymmetric models in which 
the capture and annihilation of the neutralino yields larger neutrino fluxes are 
inconsistent candidates for the primary component of the galactic halo. (To be 
precise, we do not implement the 2 GeV cutoff in our calculation, but since we are 
primarily interested in heavy neutralinos here the fraction of our signal from lower 
energy neutrinos should be insignificant.) In Fig. 9 the dark shading denotes the 
regions of parameter space excluded by this constraint. The light shaded regions 
are those that would be excluded if the observational flux limits were to be 
improved by a factor of 100. The curve inside the light shaded areas encloses 
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regions of parameter space that would be excluded if current observational limits 
were improved by a factor of 10. To indicate the sensitivity of these results 
to uncertainties in the calculation, the dashed curve inside the excluded region 
indicates the region excluded if the true neutrino rate is only 1/5 as large as our 
calculations indicate. In (a) tan/? = 2, myo =20 GeV, the squark mass is taken 
to be infinite and y > 0, and (b) is similar except that y < 0. In (c) tan (3 = 2 
and my o = 20 GeV, in (d) tan/? = 2 and myo = 35 GeV, and in (e) tan/? = 25 
and my o = 35 GeV. In (c), (d), and (e), the squark mass is assumed to be 20 
GeV greater than the neutralino mass and only regions of positive y are shown. 

From Fig. 9, we see that limits on energetic neutrino fluxes from the Sun 
already exclude many supersymmetric models with heavy mixed-state neutralinos 
lighter than about a TeV when the lightest Higgs is light and t an (3 is small 
[Fig. 9(a), (b), and (c)], or when tan/? is large [Fig. 9(e)], independent of the 
squark mass. Unfortunately, the region of m ya-tan (3 parameter space in which 
current neutrino limits might exclude neutralinos as dark matter candidates is 
similar to that excluded by current LEP results;" the rates for neutrino events 
from models with larger values of m yo [Fig. 9(d)] are much smaller. Also, current 
neutrino-flux bounds are ineffective in ruling out neutralinos that are almost pure 
Higgsino or B-ino; however, if the observational bounds are improved by a factor 
of ten, far more supersymmetric dark-matter candidates would be observable. 
For values of tan/? and my o near the current observational limits [Fig. 9(a), (b), 
and (e)], most heavy Higgsinos would be observable, independent of the squark 
mass, should they be the primary component of the galactic halo; for larger my o, 
the rates are smaller [Fig. 9(d)]. The rates from heavy B-inos are sensitive to 
the squark mass as may be seen by comparing Fig. 9(a) and Fig. 9(c). If the 
squark mass is much greater than the neutralino mass [Fig. 9(a)], then B-inos 
that are extremely pure will not be observable, but if the squark mass is near 
the neut r alino mass [Fig. 9(c)], the event rates are much greater. Also, note that 

the event rates are much smaller from supersymmetric models with negative y. 
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This is because the elastic-scattering cross sections are generally smaller " which 
leads to a smaller capture rate. 

Throughout we have taken the top-quark mass to be 120 GeV; however, our 
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results axe generally insensitive to this assumption. This is because the event 
rates are determined primarily by the capture rates in the Sun which do not 
depend on the top-quark mass. Increasing the top-quark mass would increase 
the fraction of annihilation products that are top quarks relative to the fraction 
that axe gauge or Higgs bosons, and the neutrino spectrum from top quarks 
is generally softer than that from gauge bosons. Therefore, an increase in the 
top-quark mass would result in a slightly lower event rate for models where the 
number of top-quark final states is comparable to the number of gauge-boson 
final states. 

- By comparing Fig. 9 with Fig. 6 we find that in the excluded regions the 
capture and annihilation rates are laxge enough that the number of neutralinos 
in the Sun has reached equilibrium (f© > t^). Generally, we find that current 
observational limits on energetic neutrino fluxes would have to be increased by 
about 2 orders of magnitude until neutralinos that have not yet reached their 
equilibrium in the Sun are detected 

VI. CONCLUDING REMARKS 

One of the most important questions facing particle physics and cosmology 
is the nature of the dark matter known to exist throughout the Universe and 
in our galactic halo. A well-motivated extension of the SU(3)c x SU( 2)i x 
£7(1 )y model of particle interactions is the minimal supersymmetric standard 
model. If low-energy supersymmetry exists in Nature then it is likely that the 
neutralino is the lightest supersymmetric particle. Although the neutralino was 
originally taken to be light, its mass could also lie in the 100-GeV range, and as 
unsuccessful accelerator searches push the mass scale for supersymmetry upward 
this possibility becomes more attractive. Calculations show that in much of 
parameter space the neutralino has a relic abundance suitable for solving the 
dark matter problem. Given this result, it remains to be seen experimentally 
whether neutralinos do indeed populate our halo. 

In this paper we have proposed that the presence of heavy neutralino dark 
matter be inferred through the observation of energetic neutrinos produced by 
neutralino annihilation in the Sun. Neutralinos that are primarily Higgsinos or a 
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mixed Higgsino/gaugino state are captured in the Sun by coherent elastic scat- 
tering due to light-Higgs-boson exchange off of nuclei in the Sun, and for mixed- 
state neutralinos the capture is quite efficient. If the squark is not much heavier 
than the neutralino, gauginos are captured via spin-dependent squark-exchange 
scattering off of hydrogen in the Sun. Since the masses of heavy neutralinos he 
outside the Earth’s resonance range, capture in the Earth is relatively inefficient. 

Neutralinos that have been captured in the Sun will annihilate therein and 
high-energy neutrinos will be produced by the decays of the annihilation prod- 
ucts. Calculation of the energy spectrum of neutrinos from such a source as they 
emerge from the Sun is quite involved as the cascade from the annihilation prod- 
ucts must be modeled considering, amongst other things, the effect of the solar 
medium on the shower. In addition, since the neutrinos have very high energies, 
absorption and energy loss of the neutrinos as they pass through the Sun must 
be included in the calculation. 

The most promising method of detection of these neutrinos is through obser- 
vation in underground detectors of upward-moving muons produced by the neu- 
trinos in the rock below the detector. Current limits from IMB on the number 
of such thoroughgoing muons may already be used to constrain regions of heavy- 
neutralino parameter space where the neutralino is a mixed Higgsino/gaugino 
state and with a mass less than about 300 GeV. Furthermore, in other regions 
of parameter space, where the neutralino is either slightly heavier (though still 
in the sub-TeV range) or closer to being a pure Higgsino or gaugino state, the 
predicted event rates are large enough that energetic neutrino signals may be 
observable in the near future with increased observing time or larger detectors. 
Given the enormous importance of such a discovery and the promise of obser- 
vation of such a signal from many supersymmetric dark-matter candidates, the 
search for energetic neutrinos from the Sun should be pursued. 

The final result of our calculation that was compared with experiment was the 
flux of neutrino-induced upward-moving muons; therefore, the strongest limits 
should eventually come from detectors with the largest surface area or longest 
exposure time. The current IMB 28 limits come from a detector of area roughly 
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400 ra 2 and an exposure time of about a year, and the limits from Kamiokande 

99 

II" come from a slightly smaller exposure. The next improvement should come 
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from MACRO which will have an area more than twice as large as IMB, and in 

the more distant future there may be a factor of 10 improvement in the collection 

area with a deep-sea detector. There is also the intriguing possibility of an 

increase in detector area of several orders of magnitude by looking for Cherenkov 
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radiation from energetic muons in deep antarctic ice. 

To see the prospects for discovery of dark-matter candidates via observation 

of muons induced by neutrinos from WIMP annihilation in the Sun let us consider 

the background of throughgoing muons induced by atmospheric neutrinos. The 

28 

flux of such muons (with energies larger than 2 GeV) is 

$„(£ > 2 GeV) = 0.075 m -2 yr -1 sr -1 . (28) 

Now although the angular size of the Sun in the sky is quite small and the 
detector resolution may be quite good, the angle between the muon direction and 
the direction of the parent neutrino has an intrinsic distribution with average of 
roughly 8^ ~ 15°/[£' # j/(2 GeV)] 1 / 2 , so muon tracks from within 15° of the Sim 
need to be accepted. We see that the background from an angular window of this 
size is comparable to the IMB limit of 0.0265 m~ 2 yr -1 . So additional exposure 
will improve this flux limit by providing the statistics needed to distinguish excess 
signal from background. 

Another strategy for improving the signal to noise ratio is to raise the muon- 
energy cutoff E c ^. Since the atmospheric neutrino flux decreases roughly as E~ 3 
(to be conservative) and the probability for detection of a neutrino of energy E„ 
is proportional to E„, the background event rate decreases only logarithmically 
with increasing cutoff energy; of course, this is not the whole story. Since the 
mean muon-production angle oc E^ 1 ^ 2 the size of the angular window around 
the Sun from which muon tracks must be accepted is accordingly smaller; con- 
sequently, the background event rate is proportional to (25£ ut ) -1 . On the other 
hand, most of the neutrinos from WIMPs with masses of 100-1000 GeV should 
have energies well above 10 GeV; furthermore, the detectability of energetic neu- 
trinos is proportional to the neutrino energy. So by accepting muons with energies 
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greater than 10 GeV, for example, the background is decreased by a factor of five 
while the dark-matter signal should be reduced only slightly. Of course, if such a 
cutoff is to be implemented the neutrino spectra from heavy- WIMP annihilation 
in the Sim should be more carefully determined, either through Monte Carlo or 
more detailed analytic modeling of interactions of decay products and neutrinos 
with the solar medium to determine exactly how much of the signal is lost by 
rejecting muon events with energies lower than the cutoff. 

We should mention that throughout we have assumed that neutralinos are 
the primary component of the galactic halo. Of course, if neutralinos constitute 
only a fraction of the dark matter, then the rates for detection will be lowered 
accordingly. There is also the question of whether the relic abundance of the 
LSP associated with a given supersymmetric model can is suitable to account for 
the dark matter in galactic halos. Generally, it is assumed that if the fraction of 
critical density contributed by neutralinos today is 0.025 < < 1, where h is 

the present Hubble parameter in units of 100 km s" 1 Mpc -1 , then the neutralino 
is a good dark-matter candidate. If fl^h 2 > l the relic density is too large to be 
consistent with the observed age of the Universe and if Q^h 2 ^ 0.025, the relic 
abundance is too small to make up the primary component of the galactic halo. 

Here we assume that all of the heavy neutralinos we consider are candidates 
for the primary component of the galactic halo. The relic abundance of a WIMP 
depends on its abundance in the early Universe at “freeze out” , when the an- 
nihilation rate of the WIMP falls below the expansion rate. The annihilation 
rate at any given time depends on the temperature of the Universe and the cross 
section for annihilation of the WIMP which is determined by the particle-physics 
model. On the other hand, since we have little familiarity with the conditions 
in the Universe before big bang nucleosynthesis, the expansion rate at freeze 
out cannot be reliably predicted. If one makes the simplest — and standard — 

assumption, that the early Universe was radiation dominated, then it is found 
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that the relic abundance of heavy neutralinos is generally greater than 0.001. 

However, many nonstandard scenarios accomodate an expansion rate at freeze 

out larger than that in the radiation-dominated Universe, so if the standard 

calculations find a relic abundance greater than 0.001, nonstandard scenarios 
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allow for a relic abundance greater than 0.025. Conversely, if standard calcu- 
lations yield Q^h 2 > 1, a value of ft*/* 2 < 1 is possible if the abundance was 
diluted by some entropy-producing process such as inflation, a quaxk/hadron or 
electroweak phase transition, or out-of-equilibrium decay of a massive particle. 
Therefore, since the standard calculations yield relic abundances for LSPs within 
a few orders of magnitude of the dark matter window, 0.025 ^ ^ 1, and the 

abundance of a thermal relic in nonstandard cosmological models may differ from 
that in the standard radiation-dominated Universe by a few orders of magnitude, 
almost all heavy neutralinos should be considered dark-matter candidates. 

Given that energetic neutrinos from heavy neutralino annihilation in the Sun 
may be observable, we speculate that neutrinos from annihilation of other heavy 
dark-matter candidates (such as Majorana neutrinos) may also be observable. 
Such a heavy WIMP would have to be captured readily in the Sun, either by 
a coherent interaction with heavy nuclei or by a sizable spin-dependent elastic 
scattering cross section that could result from the exchange of another particle 
not much heavier than the WIMP (t.g., a heavy lepton in the case of a Majorana 
neutrino), or maybe by a strong coupling to the Z. Even if the dark matter con- 
sists of some heavy WIMP other than the MSSM neutralino, the MSSM provides 
a good example of a particle-physics model with a well-determined phenomenol- 
ogy that is consistent with current laboratory results and contains am excellent 
dark-matter candidate. This example shows that the idea that galactic halos are 
populated by (possibly detectable) WIMPs is alive and well and that the quest 
for their discovery should be pursued vigorously. 

To conclude, we note that the properties of the heavy neutralino in many 
models are such that their capture and annihilation in the Sun yields an observ- 
able flux of energetic neutrinos. We also point out that in many models, a heavy 
neutralino may easily make up the primary component of the galactic halo while 
remaining invisible to neutrino detectors, so null results from energetic neutrino 
searches are not likely to rule out supersymmetric dark matter. Nevertheless, 
given the present uncertainty as to the nature of the dark matter, the popularity 
of supersymmetry in particle physics, and the interesting “coincidence” that the 
relic abundance of the LSP in most supersymmetric models falls near the dark- 
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matter window, it is clear that the search for energetic neutrinos from the Sun 
holds considerable promise for discovery, should neutralinos reside in the galactic 
halo. 
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APPENDIX A: ELASTIC SCATTERING CROSS SECTION 

The neutralino may elastically scatter off of a nucleus via a scalar interaction 
where the WIMP interacts coherently with the entire nucleus, and if the nucleus 
has spin the neutralino may also scatter via an axial interaction. The cross 
section for scattering of a neutralino off of nucleus t via an axial interaction (the 
“spin-dependent cross section”) is 5 

2 


<7SD = 




7 r(m^ + m,) 2 3 




I u t d,s 


where 


4 = - z&) 

2mlcP„ 


- x 2 ^ + [7* L Z i2 ~ tan 6 W (T* L - e,)Z tl ] 2 + tan 2 d w e 2 q Z 2 x \ , 

(A2) 


l w 


and 


0 

X" = 


m 


w 


q (m j + m + i) 2 — (Mq — m,) 2 1 


(A3) 


36 


(A4) 


31 


is the squark-exchange suppression factor, and 


A = i {1 + + !)-<(»+ Dl/W + 1)]} . 

is the Lande factor from the one-particle nuclear shell model for a nucleus with 
spin J and an unpaired nucleon with spin s p and orbital angular momentum l. 
Here, m q is the (current) quark mass, d q = — Z l3 / cos/? for down-type quarks, 
d q = Zi 4 /sin/3 for up-type quarks, T% L is the weak isospin of the quark, e q is its 
charge, and $w is the Weinberg angle. The quantity A q measures the fraction 
of the nucleon spin carried by the quark. In the naive flavor-SU(3) quark model 
Au = 0.97, Ad = —0.28 and As = 0; however, the EMC collaboration reports 
Atz = 0.746, Ad = -0.508 and As = -0.226. 36 

For the capture-rate calculation the spin- dependent cross section and the sum 
that appears in Eq. (Al) may be simplified considerably. The only element with 
spin in the Sun found in abundance is hydrogen. For hydrogen (4/3)A 2 J( J + 1) = 
1 find in the EMC model 


E — 0.37(Z, 2 3 — Zj) — x 2 


Z 2 Z 2 

- 3.98 x 10 -7 — £-5 + 2.86 x 10 _9 ^f- 
cos 2 ft sm" ft 


+ 0.003Z 2 2 + 0.133Z, 2 Za -1- 0.073^ , 

(A5) 


while in the flavor-SU(3) model 
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cos 2 ft 
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The term proportional to (Z 2 3 — Z 2 4 ) arises from Z exchange, and the second 
term arises from squark exchange. For heavy 5-inos Z { 3 ~ Z,4 ~ 0, for heavy 
Higgsinos Z 2 3 ~ Z 2 4 , and as we will see below, for heavy mixed-state neutralinos 
the axial interaction is much weaker than the coherent interaction; therefore, 
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scattering of heavy neutralinos via Z exchange is essentially negligible. In addi- 
tion, from Eqs. (A5) and (A6) one can see that if the neutralino is pure Higgsino 
spin-dependent scattering due to squark exchange is also negligible, but if the 
neutralino is pure B - ino (Z , u ~ 1 and Zn ~ Z,3 ~ Z,4 ~ 0) and the squark is not 
much heavier than the neutralino then spin-dependent scattering due to squaxk 
exchange may be significant. By comparing Eqs. (A5) and (A6) we also see that 
had we used the flavor-SU(3) quark model the capture rates would be roughly 
3 times as large as those obtained using the EMC results which we used in this 
work. 


The cross section for scattering via a scalar interaction is obtained from 
Refs. 5, 32, 33, and 21. Griest 5 obtained the results for a coherent scalar in- 

32 

teraction via exchange of a virtual squark, and Barbieri, Frigeni, and Giudice 
obtained results for a coherent scalar interaction in which a Higgs boson is ex- 
changed; however, in both of these papers it was assumed that the nucleon mass 

is due to gluons. 48 Recent measurements of the pion- nucleon sigma term imply 
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that a significant fraction of the nuclear mass is due to a sea of strange quarks. 
When applied to coherent neutralino-nucleus scattering it is found that although 
the component of squark- and Higgs-nucleon coupling due to gluons is reduced, 
there is an additional component due to squark and Higgs coupling to the strange- 
quark sea and the net effect is a significant increase in the squark- and Higgs- 
nucleon coupling. 21 


The scalar cross section may be derived from the effective Lagrangian 


5,32 
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where k [ ^ = sin a/ sin/? and k = cos a/ sin for up-type quarks, = 
cos a/ cos /? and k^ = — sin a/ cos 0 for down-type quarks, gu 7 = (Z,3 sin a + 
Zt { 4 cos 0), and gn x = {Zi^ cos ot + sin a). In addition to terms due to exchange 
of the lightest Higgs boson and the squark, to be complete we have included 
a term due to exchange of the heaviest Higgs boson although it should generally 
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be smaller than that due to exchange of the lightest Higgs boson. 

The scattering cross section is obtained from the square of the matrix element 
(f\C e fi\i) of this effective Lagrangian between the initial and final neutralino- 
nuclear states. In Ref. 48 (as modified by Ref. 33) it is shown that the coupling 
of a scalar field to the gluons in the nucleus occurs via a heavy-quark (c, b, and 
t quarks) loop so that 

(N\m k hh\N) = ^m,(0.56), (A8) 

where h is a heavy-quark field, m/, is the heavy-quark mass, and \N) is the nuclear 
wave function. In addition, measurements of the pion-nucleon sigma term imply 
that 33 

(N\m s Ss\N) = ^m t (5.94), (A9) 

where s is the strange-quark field and m s is its mass. The matrix elements of 

m q qq for the u and d quarks are much smaller. With these results it is easy to 

find that the matrix element is 
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and the cross section for scattering off of nucleus i via a coherent scalar interaction 
is 


<7SC = 


4m * 2m L|</|£ effl o p. 


nfax + m t ) 2 


(All) 


We should clarify that this is the cross section that would be measured only 
if the neutralino interacted coherently with the entire nucleus. If the inverse 
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of the momentum transfer 1 /q in the scattering event is small compared with 
the nuclear radius R then the neutralino does not interact coherently with the 
entire nucleus and the actual cross section is momentum-transfer dependent (or 
equivalently, scattering-angle dependent) and is given by Eq. (All) times |F(g 2 )|, 
the form-factor suppression. The effect of the form-factor suppression on capture 
in the Sun and Earth is discussed in Section II. 

APPENDIX B: MIXED GAUGE/HIGGS BOSON FINAL STATES 

In addition to the gauge-boson and Higgs-boson final states considered in 
Ref. 6, neutralinos may annihilate into mixed Higgs/gauge boson final states 
when the mass of the neutralino exceeds the average of the gauge- and Higgs- 
boson masses. At zero relative velocity the available channels are ZH®, ZH ®, 
W + H~ , and W~ H + . Annihilation into ZH ® is possible in general, but does not 
occur at zero relative velocity for CP-conserving theories. The reason is that 
at zero relative velocity, neutralino-neutralino annihilation occurs via an 5 -wave 
and due to Fermi statistics, the initial state has CP = —1. Since the Z has spin 
1 and the Higgs is a scalar, the orbital wave function of the outgoing state must 
have / = 1, and since the Z is CP-even and the H 3 is CP-odd, the final state 
must have CP = 1 and is therefore inaccessible from the initial state. 

Since the ZH® final state is the first mixed channel to open up as the neu- 
tralino mass is increased, we will consider it first. [Incidentally, since (m# 0 + 
mz)l 2 may be less than m^, this channel may be open for neutralinos that 
are lighter than the W, a possibility that was not considered in previous work.] 
Throughout this Appendix we will use the notation of Griest, Kamionkowski, and 
Turner (GKT)f and some of the couplings we will use here are defined there. 

Annihilation of two neutralinos into ZH ® occurs via s-channel exchange of a 
Z and a and by t- and u-channel exchange of all four neutralinos. The cross 
section crzH° f° r this process as relative velocity v re i — *■ 0 is 


a ZH$ v rel 


327rm^ 3 ’ 


(Bl) 


40 



where 
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is the momentum of the outgoing paxticles and 
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Here, z = mz sin(/3 — a)/ cos Q\y is the coupling at the H^ZZ vertex, F,y = 
(2,3 Z j3 - ZhZj4)/2cos0w is the coupling at the Zx®x® vertex, A/,-,* is the 
coupling and is given in Eq. (C9) of GKT, h = cos(a — /3)/2cos Bw 
is the ZH^H® coupling, the sum is over all four neutralinos, and t = [( m | + 

m H° )/2] — m x 2 - 

For larger neutralino masses the Z channel opens up. (Recall that the H j 
is always heavier than the 2.) The cross section for annihilation into ZH J* may be 
obtained from that for annihilation into ZiJ® by simply replacing mg o by mg o, 
M 2 ij by M\ij, and using z = mz cos(/3 — a)/2 cos^ra and h = sin(ar — /?)/2cos#^. 

Annihilation into WH ± final states occurs through s- channel exchange of the 
H® and t- and u-channel exchange of the two charginos. The cross section for 
this process as relative velocity t? re i — ► 0 is 
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The sum is over the two charginos, and the quantities e, and f x are given in GKT 
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Eq. (A2), and 
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where the angles <j>+ and <j>- are related to the diagonalization of the chargino 
mass matrix and axe given in Ref. 49, and e = det X /\ det X| and X is the matrix 
defined in Eq. (C9) of Ref. 3. Here, 
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and t = [( m^r + m# + )/ 2] — m * 2 . 
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FIGURE CAPTIONS 


1. Lightest neutralino composition and mass for tan /3 — 2. The broken curves 
are contours of constant neutralino mass and the solid curves are con- 
tours of constant gaugino fraction ( Z 2 j 4- Z 2 2 ); in (a) /i > 0 and in (b) 
fi < 0 . 

2. Form-factor suppression of the rate of accretion of heavy WIMPs onto the 
Sun from scattering off of nuclei with atomic masses 4, 12, 16, 24, 32, and 
56 as a function of the neutralino mass, 

3. Contour plots of the capture rate of neutralinos in the Sun assuming neu- 
tralinos make up the primary component of the dark matter and that the 
squark mass is infinite. The double curve indicates a capture rate of 10 24 
s -1 ; the spacing between other curves are decades, the capture rate decreas- 
ing toward higher masses. In (a) tan/3 = 2, mjyo =20 GeV, and n > 0 
and (b) is the same except /a < 0. In (c) tan/3 = 2 and m#o = 35, and in 
(d) tan /3 = 25 and m H o = 35. In (c) and (d) only regions of positive (i are 
shown; the plots for negative /i are similar. For convenience, the mass and 
composition contours are also shown. 

4. Same as Fig. 3(a) but here the squark mass is assumed to be 20 GeV heavier 
than the neutralino mass. 

5. Contours of the fraction of the capture rate due to spin-dependent scatter- 
ing when the squark is assumed to be 20 GeV heavier than the neutralino, 
and tan/3 = 2 and m H o = 20. In the shaded regions the fraction is greater 
than 0.5, and the contours indicate where the fraction is 0.01, 0.5, and 0.99. 
Again, mass and composition contours are also shown, and plots for other 
values of tan/3 and o are qualitatively similar. 

6. Contours of t©/r a- In the dark shaded regions Iq/ta < 0.33 and in the light 
shaded region t©/r a < 1-82; elsewhere, t©/^ > 1.82. In (a) tan/3 = 2, 
mu o = 20, the squark mass is taken to be infinite, and fi > 0; (b) is similar 
but /x < 0 is shown; and (c) is similar to (a) but the squark mass is taken 
to be 20 GeV heavier than the neutralino mass. Plots for other values of 
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tan/3 and m# o axe similar. 

7. Second moments of the neutrino yields from the Sun from the cc, bb, tt, 
r ± , W ± , and H°H° (where tan /3 = 2 and m H o = 20 GeV) annihilation 
channels as a function of the neutralino mass. 

8. Contours of the fraction of the neutrino signal that comes from gauge-boson 
final states. In the shaded regions the fraction is greater than 0.5, and the 
contours indicate where the fraction is 0.01, 0.5, and 0.99. In (a) the squark 
mass if taken to be infinite, and in (b) the squark mass is assumed to be 20 
GeV heavier than the neutralino mass. In both, tan/3 = 2 and m#o = 20 
GeV and fi > 0. Plots for other values of tan/3 and m# o and for negative 
y. are qualitatively similar. 

9. Regions where the neutralino is excluded as the primary component of 
the galactic halo by limits on the flux of upward-moving neutrino-induced 
muons from the Sun. The dark shaded regions are those excluded by current 
IMB limits. The light shaded regions axe those that would be excluded if 
current observational limits were improved by a factor of 100. The curve 
inside the excluded region encloses the region that would be excluded if the 
true neutrino flux was 1/5 of the results of the calculation here, and the 
curve inside the light shaded region encloses regions that would be excluded 
if the current observational limits were improved by a factor of 10. In (a) 
tan /3 = 2, m H o = 20 GeV, n > 0, and the squark mass is taken to be 
infinite and (b) is the same except /i < 0. In (c) tan/3 = 2 and o = 20 
GeV, in (d) tan/3 = 2 and m H o = 35 GeV, and in (e) tan/3 = 25 and 
m H o = 35 GeV. In (c), (d), and (e), the squark mass is assumed to be 20 
GeV greater than the neutralino mass and only regions of positive n axe 
shown. Plots for negative n axe similar, but excluded regions are smaller. 
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